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ABSTRACT 



A channel with rectangular cross-section, 40 to 1 aspect ratio (height is 0.0127 in) 
and 4.27 m test section length (336 channel lengths) is used to study the elfects of im- 
posed pulsations on laminar, transitional, and turbulent How phenomena. Periodic \e- 
locity variations arc produced in the test section using a single rotating vane located in 
the How downstream of the test section. Survey data for ^’/d of 0.50, 0.85, and 0.90 
show how normalized average velocity and turbulence intensity data vary with Reynolds 
number for Reynolds numbers ranging from 1 100 to 3400 at Stokes numbers of 4.08 and 
5.79 and Strouhal numbers from 0.033 to 0.121. These conditions are produced from 
pulsations imposed at 1 Ilz and 2 llz to show how transition develops both with and 
without imposed pulsations. Intermittcncy variations increase from near zero to one as 
Reynolds number changes from 1450 to 2100 for Y/d = 0.90, 0.85, and 0.50. Time- 
averaged magnitudes of the normalized longitudinal turbulence intensity increase signif- 
icantly as Re changes from 1450 to 3100 at Y/d = 0.90. A maxima is reached near 
Re= 2400. Afterwards these data show a decrease and level off in the turbulent region. 
At Y/d = 0.85, the same trend in normalized longitudinal turbulence intensity appears 
with a maxima appearing near Re= 2300. Tor the center of the channel (Y/d = 0.50), the 
maxima appears near Re =2000. The normalized longitudinal turbulence intensity level 
in the turbulent region increases with Y/d. Of particular interest are conditions where 
the normalized longitudinal turbulence intensity is greater near the center of the channel 
than near the walls. This occurs for Reynolds numbers of 1500 to 2500. Spectra show- 
evidence of intermittency and drastic increases in fluctuating energy starting at Re= 1500 
and continuing to Re= 1700. 
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I. INTRODUCTION 



A. BACKGROUND 

Flow pulsations and unsteadiness are present in many practical engineering flow 
situations. Most flows of technological interest also undergo transition from a laminar 
to turbulent state. Transition is by itself extremely complex; however, with flow 
pulsations and unsteadiness, accurate prediction of the location and extent of transition 
and accompanying changes of important flow properties is even farther beyond present 
computational modeling abilities. Consequently, experiments are necessary- to elucidate 
flow behavior. According to Shemer [Ref 1], such experiments should be "restricted to 
the simplest and most well-defined conditions" to "obtain a clear physical description of 
a fluid dynamical process." 

The present study employs a straight channel to investigate a situation consistent 
with Shemer's suggestion. This flow is also studied because little experimental informa- 
tion is available for transitioning flow in straight channels with imposed pulsations. 
Some experimental results are published for pipe flows as well as limited analytic and 
numeric results for channel flows; however, the author knows of no other experiments, 
other than the continuing current work, performed to study the stability or transition 
of plane channel flow with imposed pulsations. 

Here, pulsations are induced with a single rotating vane located downstream of the 
flow test section. The device is similar to the rotating profiled sleeves employed by Tu 
and Ramaprian [Ref 2], Ramaprian and Tu [Ref 3] and Stettler and Hussain [Ref 4] to 
study the influences of sinusoidal pulsations on pipe flows. Their sleeves are used in 
water and operate simply by altering the exit flow area. Other studies of influences of 
imposed pulsations on pipe flows use reciprocating pistons in cylinders to produce 
pulsations [Refs. 5,6,7, 8 ,9,10, 11 ]. Simpson et al [Ref 12] use a rotating blade damper 
system at the inlet of their wind tunnel, upstream of screens and honeycomb, in order 
to study the effects of pulsations on turbulent boundaiy layer flows. By changing the 
motion of the five individual damper blades, sinusoidal perturbations to the flow are 
produced with amplitudes amounting to 11-93 % of the maximum velocity over fre- 
quencies from 0.6-2 Hz. 

Of the experimental studies of transitional flows in pipes, .Merkli and Thomann [Ref 
5], Sergeev [Ref 6], Hino and Sawamoto [Ref 7], Hino et al [Ref 8], Gerrard [Ref 9], 
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Gilbrcch and Combs (Ref. 10), Sarpkaya [Ref. 11], Ramaprian and Tu (Ref. 3], Shemer 
[Ref 1] and Stettler and Hussain [Ref 4] also examine the effects of imposed pulsations. 
Of these studies, [Refs. 5-8] investigate the instability of sinusoidally modulated pipe flow 
with zero mean flow. Gerrard [Ref 9] conducts his experiments at transition Reynolds 
number of 3770. He relates variations of turbulent flow over individual flow pulsations 
to magnitudes of vorticity diffusion. Gerrard also suggests that, with high pulsation 
frequency, viscous effects are confined to a very thin Stokes layer near the wall where 
the fluid is retarded. Gilbrech and Combs [Ref 10] and Sarpkaya [Ref 11] examine the 
effects of amplitude and frequency of imposed pulsations on the gro\nh rate of arti- 
ficially introduced plugs. Both investigators indicate that pulsations increase the critical 
Reynolds number as long as they are not so large as to cause local flow reversals. 
Ramaprian and Tu [Ref 3] study transitional pipe flow at a mean Reynolds number of 
2100. They indicate that flow pulsations at frequencies ranging from 0.05 to 1.75 Hz 
increase the critical Reynolds number only when the turbulent intermittency is small. 
With higher intermittency levels, the authors suggest that imposed pulsations affect the 
flow only when the pulsation frequency is near the characteristic frequency of the tur- 
bulence. Shemer's [Ref 1] work focuses on large pulsation amplitudes at a single low 
frequency. For a mean Reynolds number of about 4000 and a pulsation frequency of 
0.37 Hz, he concludes that transition is governed primarily by the instantaneous 
Reynolds number. Stettler and Hussain [Ref 4] present results for Stokes numbers as 
high as 70. The authors provide a three-dimensional map of stability-transition regimes, 
and also indicate that transition in pipes is associated with plugs of turbulence which can 
grow or shrink in size. 

Numerical and analytical investigations of the influences of imposed pulsations on 
flow in a plane channel are described by Grosch and Salwen [Ref 13], Herbert [Ref 
14], Hall [Ref 15], von Kerezek [Ref 16] and Singer, et al [Refs. 17,18]. In addition to 
these studies, Tozzi [Ref 19] provides results of a numerical investigation of the influ- 
ences of imposed pulsations on flow in a pipe, and Davies [Ref 20] provides a review 
of much work done on the stability of time periodic flows. Of the above mentioned 
analytic studies of plane channel flow, Grosch and Salwen [Ref 13] solve a set of 
linearized equations by integrating through one period of pulsation. They conclude that 
flows with small amplitudes of imposed pulsations arc more stable than steady flows, 
where the degree of stabilization depends upon interactions between shear waves gener- 
ated by the imposed pulsation and flow disturbances. A modified version of Grosch and 
Salwen’s energy analysis is described by Herbert [Ref 14]. Hall [Ref 15] reports results 
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obtained with high frequency pulsations and concludes that such pulsations are slightly 
destabilizing irrespective of amplitude, von Kerczek [Ref 16] modulates the pressure 
gradient in a perturbation analysis, and shows that the oscillating flow is more stable 
than the steady flow for frequencies of imposed pulsations greater than about one tenth 
of the frequency of the steady flow neutral disturbance. However, at both very low and 
very high frequencies of imposed pulsations, the flow is slightly unstable. Singer et al 
(Refs. 17,18] describe results from a direct Xavier-Stokes simulation of flow in a plane 
channel with imposed pulsations and indicate that imposed sinusoidal pulsations provide 
a stabilizing effect at all but ver}’ low frequencies. Significant variations in the ampli- 
tudes of Tollmein-Schlichting waves are also noted which depend upon the Strouhal 
number and the amplitude of the pulsations, as well as on initial amplitudes of the 
Tollmein-Schlichting waves. Tozzi's (Ref 19] study of pipe flo\v shows that imposed 
pulsations are stabilizing up to very high amplitudes of imposed pulsation. 

B. OBJECTIVES 

This study provides experimental data which show the influences of imposed 
pulsations on laminar, transitional, and turbulent flows in a rectangular cross-section 
channel with 40 to 1 aspect ratio and L/d of 336 (d = 0.0127m). The imposed pulsation 
device employed is based on existing designs, but here is applied to channel flow with 
air. Mean streamwise velocity and longitudinal turbulent intensity velocity data as a 
function of Reynolds number for various pulsation frequencies are given for three lo- 
cations across the channel height. The current study is complementary to the wo-^’k of 
Morrow (Ref 21] and is an extension of work by Koth (Ref 22], Coumes (Ref 23], and 
Longest (Ref 24]. These previous studies found imposed pulsations at 1 and 2 Hz to 
have a greater destabilizing effect on the longitudinal velocity fluctuations in the Stokes 
layer at some Reynolds numbers than observed in the present study. The present results 
are given for Strouhal numbers from 0.039 to 0.121, and Reynolds numbers ranging from 
1100 to 3400 (Vjrt from 1.3 to 4.1 m,'s), where each is based on bulk mean velocity and 
channel height. To produce these conditions, pulsations are imposed on the flow at 
frequencies of 1 Hz and 2 Hz. These imposed pulsation frequencies correspond to 
kinematic viscosity based Stokes numbers of 4.08 and 5.79, respectively. 
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C. THESIS ORGANIZATION 

In the chapters that follow, the experimental facilities, including the channel, un- 
steadiness device and data acquisition system, are detailed in Chapter II. The exper- 
imental procedures for calibration of the hot-wire probe, channel validation, flow 
measurement and spectral analysis are described in Chapter III. Results are discussed 
in Chapter IV. A summarx’ of results and conclusions appear in Chapter V. 
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II. EXPERIMENTAL FACILITIES 



A. CHANNEL 

The straight channel was specifically designed for the present study and is located 
in the laboratories of the Department of Mechanical Engineering of the Naval Post- 
graduate School. Figure 1 is a schematic of the channel, and the coordinate system 
employed, x is the longitudinal coordinate measured from the downstream edge of the 
nozzle, y is the normal coordinate measured from the bottom wall, and z is the trans- 
verse coordinate, measured from the spanwise centerline of the channel. 

The test section of the channel is 4.27 m long, 1.27 cm high and 50.8 cm wide (inside 
dimensions). This gives an internal aspect ratio of 40 to 1. The top and bottom walls 
are 6.35 mm thick lexan supported by ribs and cross beams along their length. The side 
walls are removable for access to the inside of the channel. Inlet flow is managed with 
a honeycomb, screens and a 20:1 contraction ratio nozzle. Two layers of cheese cloth 
are also placed at the inlet to remove dust and dirt from the entering air. At the test 
section exit are three 10.16 cm long frames. Screens are between the flanges of the 
frames and honeycomb is placed just upstream of the last screen. The unsteadiness de- 
vice is housed in the middle frame. A two-dimensional diffuser (45.72 cm long, 3 degree 
total angle), located downstream of the honeycomb section, connects the straight section 
to the exit plenum. This plenum provides a uniform volume of low pressure air at the 
diffuser exit. The plenum is a cube with inside dimensions of 60.96 cm along each side 
and is connected to the suction side of a 5 hp blower via piping of 5.08 cm inside diam- 
eter. Flow through the channel is metered by a 3.81 cm orifice plate assembly located 
near the plenum exit. The channel is rigidly mounted on steel frames with supporting 
ribs along its length. With this mounting arrangement, no measurable deflection of the 
channel walls occurs with either steady or pulsating flow. Vibration damping systems 
are employed on the mounting frame and test section to further reduce vibrations due 
to mechanical sources. Additionally, all re.sults are obtained at night with no other ac- 
tivity in the laboratory to minimize disturbances from acoustic sources. 

Channel bulk flow velocities ( F*,*) of 15 m,'s are possible; however, the present study 
uses a range of F*,* from 1.3 to 4.1 nrs. The channel is designed so transition will occur 
after laminar flow has developed fully. The measured turbulence intensity of the flow 
at the nozzle exit (Figure 7) is 0.12 to 0.14 percent for Reynolds numbers less than 2100 
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with no imposed pulsations, as well as with pulsations imposed at 1, 2, 3, and 4 Hz. 
As Re increases to 3600, inlet intensities both with and without imposed pulsations are 
slightly higher (0.17 to 0.24 percent). 

B. UNSTEADINESS DEVICE 

Periodic variations of the flow rate can be induced in the test sections of the wind 
tunnels, channel flows and pipe flows by introducing a perodic blockage at the inlet or 
exit of flow passages. The most important requirement of such an unsteady device is the 
ability to produce controlled deterministic and periodic unsteadiness without adding 
other disturbances. In an open circuit induction channel, the unsteady device is best 
located just downstream of the test section, thereby eliminating convection of wakes and 
other flow disturbances into the test section. 

A single rotating vane is used to impose pulsations in the present study, as shown 
in Figure 2. A single vane makes this possible since the interior height of the channel 
is relatively small (1.27 cm). The vane is a 3.2 mm thick brass strip with rounded edges, 
and is supported at the ends by a 3.2 mm diameter shaft with bushings fitted to the side 
walls of the frame. Three intermediate spanwise struts increase the rigidity of the vane 
as it rotates. One of the shafts extends to accommodate a 48 tooth spur gear driven by 
a 12 tooth spur pinion mounted on the shaft of a Superior Electric, M092-FD310 
Stepper Motor. The motor is driven by a Modulynx MITAS PMS085-D050 Drive 
controlled by a Modulynx MITAS PMS085-C2AR Drive Controller. This arrangement 
allows the motor shaft to be positioned at l,'200th increments of one motor shaft revo- 
lution. The 1 to 4 gearing ratio, coupled with two cycles of imposed flow pulsations per 
vane rotation, gives one motor revolution corresponding to one half of one cycle of an 
imposed flow pulsation. The imposed pulsation frequency is set by programming the 
Drive Controller for a particular motor speed which then sets the vane rotation rate. 
The amplitude of the imposed pulsations is altered by changing the width of the vane. 
A vane width of 8.7 mm is employed in the present study. 

When the vane rotation rate is small enough, the vane produces a quasi-steady flow 
and the flow behaves as if the vane were fixed at each position and allowed to reach 
steady state behavior. The variation of the mean velocity under such circumstances is 
approximately sinusoidal, as illustrated in Figure 3. This waveform results because, 
momentarily, mass flow rates (and resulting velocities in the straight test section) are 
roughly proportional to the flow passage area provided by the vane as it rotates. In 
Figure 3, flow resistance is ma.ximum and velocity in the straight section is minimum 
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when the %ane is normal to the flow with position 2. With vane position 4, the flow 
resistance is minimum and the velocity is maximum. 

Actual flow around the rotating vane is quite complex. Velocity and pressure vari- 
ations are affected not only by the varying flow resistance from the vane, but also by 
other fluid dynamic effects such as dynamic flow separation and flow inertia. At large 
vane speeds, the vane rotations will likely result in a local swirl with blockage like a solid 
cylinder and minimal bulk flow periodicity. Imposed oscillatory' flow behavior is thus 
different and difficult to predict compared to quasi-steady flow. The assessment of the 
performance of the rotating vane as it produces imposed velocity pulsations in the test 
section is given by Ligrani, et al [Ref 25). 

C. FLOW MEASUREMENT INSTRUMENTS 

Figure 4 presents a schematic of the flow measurement equipment. 

1. Hot-«ire Anemometer 

A DANTEC 55P51 hot-wire probe, with sensor diameter of 5/j.m and sensor 
length of 1.25 mm, is used for instantaneous velocity measurements. The probe is 
mounted with the wire horizontal and normal to the flow direction. The probe position 
is controlled by a rotatable lever arm which can adjust through a range of Y,'d of 0.05 
to 0.95. The lever arm is calibrated for probe positioning with an accuracy of approxi- 
mately 0.5 mm with respect to the y coordinate. For all measurements presented, the 
probe is located 3.74 m from beginning of the test section and 97.6 mm left of spanwise 
center. These are equivalent to x,'d = 294.5 and z,'d= 7.7. 

2. Hot-Mire Bridge, Signal Conditioner, High Speed Data Acquisition System and 
Data Storage 

A DISA 55.\110 constant temperature bridge is used to operate the hot-wire 
probe at an overheat ratio of 1.8. When connected to this bridge with the same settings 
used for measurements, the hot-M'ire is calibrated in the potential flow of a wind tunnel. 
During measurement, the DC voltage from hot-w'ire bridge is measured using a 
Hewlett-Packard 3466A digital multimeter. 

A DANTEC Model 56N20 signal conditioner is used to amplify and filter the 
voltage from the bridge. During measurements, an amplifier gain of 10 is used, the 
low-pass filter is set to lOkHz and the high pass filter is set at 0.1 Hz to remove high 
frequency noise signal as well as the DC signal. The output of the signal conditioner is 
fed to a Hewlett- Packard 6944 Multiprogrammer with a buffered high speed A,'D con- 
version system. The .Multiprogrammer cards are driven using Hewlett-Packard 
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CAT 14752A software on a Hewlett-Packard 9000 scries Model 310 computer. 
Digitized instantaneous voltages are stored on 3.5 inch diskettes for further processing. 
The heading to each data set includes the time, date, bulk velocity, pulsation frequency, 
and average voltage from the digital multimeter. 

3. Phase-Averaging 

The motor used to drive the vane generates voltage pulses sequentially in each 
of the four field windings used to step the motor. Fifty pulses are available to each field 
winding for each revolution of the motor shaft, giving a total of 200 pulses per revo- 
lution. The pulses from one of the four field windings are used to trigger the data ac- 
quisition and are precisely synchronized with the vane rotation. Thus, for a given vane 
imposed pulsation frequency, f„^, the data sampling rate (TR) is given by: 

TR(samples. sec) = 

f„,, X (pulses, 'motor shaft rev.) x (motor shaft rev./imposed pulsations) 

= f„,. X 50 X 2 
= 100 X 

With this approach, 100 samples of data are obtained, spaced uniformly, over each pe- 
riod of flow pulsation, providing a phase resolution of I/lOO of the imposed pulsation 
cycle. 

Data are acquired for 640 sequential How pulsations so that the memory buffer 
of the multiprogrammer (64 kbytes) is completely filled. Afterwards, data are stored in 
the computer memorv', where they are packed into 640 arrays of 100 samples each. Data 
are then stored on a disk for later processing. Because non-uniform width voltage pulses 
may be generated by the motor drive as the motor is started and stopped, the first 6000 
points are ignored as data are processed, and all acquisitions are completed well before 
the vane motor is stopped. When no pulsations arc imposed on the flow, data acquisi- 
tion is triggered using a signal from a square wave function generator set to the same 
frequency as employed with imposed pulsations. 
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III. EXPERIMENTAL PROCEDURES 



A. HOT-WIRE CALIBRATION 

The hot-wire is calibrated in the freestream of the wind tunnel located in the labo- 
ratories of the Mechanical Engineering Department of the Naval Postgraduate School. 
During calibration, the hot-wire probe is mounted normal to the flow in potential 
freestream flow and connected to the hot-wire bridge discussed previously. Output from 
the hot-wire bridge is read on a Keithly 169 Digital Multimeter. Freestream velocity in 
the wind tunnel is measured utilizing a Kiel pressure probe, a wall static pressure tap, 
and a Validyne PS 309 Digital Manometer. Voltage and differential pressure readings 
are taken for a range of pressure drops corresponding to velocities between 1.0 and 4.0 
m,'s. Following the calibration, the hot-wire is installed in the straight channel. 

B. DATA PROCESSING 

1. Bulk Velocity 

The bulk flow velocity is determined from the pressure drop across the orifice 
plate located in the outlet piping which connects the outlet plenum and the blower. The 
channel Reynolds number and channel Strouhal number are then determined using this 
bulk flow velocity. The pressure drop is measured using a Validyne Model PS 309 Dig- 
ital Manometer. Figure 5 shows the relationship of bulk velocity to pressure drop. 

2. Instantaneous Velocity and Phase Averaging 

Instantaneous voltages are converted into instantaneous velocities using a 
look-up table that accounts for hot-wire calibration coefficients, amplifier gains and 
mean voltage levels. WMth imposed periodic flow, instantaneous velocities are given by 
U = u-fu-fu' , where u is the time-averaged velocity, u is the periodic velocity, and u' 
is the fluctuating component. 

In the present study, u and u are combined as the phase-averaged velocity 
(u), such that U = u-fu', following Ramaprian and Tu (Ref. 3]. Thus for a steady flow 
with no periodic velocity, u is then equal to u. u is then determined from phase- 
averaging instantaneous velocity results using the equation given by : 

1 ncycl 
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Here, m and n correspond, respectively, to the number of cycles and locations across 
each phase where data are sampled. Typically 5S0 cycles are ensemble averaged to ob- 
tain a phase-averaged velocity trace. 

C. SPECTRAL MEASUREMENTS 

Spectra of analog voltages from the hot-wire bridge are determined using a 
Hewlett-Packard Model 3562A Dynamic Signal Analyzer (DSA). Spectral coordinates 
are log ((volts)Vfrequency (Hz) versus frequency (Hz)), where frequency varies from 0 
to 200 Hz. Typical settings of other parameters are shown in Figure 31. With this ar- 
rangement, five sequential spectra are ensemble averaged to produce final results ob- 
tained for Reynolds numbers between 1400 and 1700. 
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IV. DISCUSSION OF RESULTS 



A. CHANNEL FLOW DEVELOPMENT LENGTHS 

Streamwise pressure distributions are shown in Figure 6 for Re= 1100, 2010, and 
3600. Streamwise pressure gradients are constant and invariant with streamwise distance 
for x'd > 55, x'd > 105, and x'd > 170 for Reynolds numbers of 1 100, 2010, and 3600, 
respectively. Such behavior evidences fully developed flow, which, if present, gives nor- 
malized velocity profiles which are self-similar with respect to streamwise distance. 
These results are in close agreement with empirical values for duct flow given by Han 
(Ref 26]. Schlichting (Ref 27], using laminar theory, gives x/d development lengths of 
44 for Re= 1 100 and 146 for Re= 3600. 

B. REYNOLDS NUMBER SURVEY OF IMPOSED PULSATION EFFECTS 

Reynolds number surveys of mean streamwise velocity and longitudinal velocity 
fluctuations are presented for Wd values of 0.50, 0.85, and 0.90 both with and without 
imposed pulsations. These surveys are based on the same profile data presented by 
Morrow (Ref 21] at x'd = 294.5 and z;d= 7.7. 

1. Local Mean Streamwise Velocity 

Figure 8 through Figure 10 show' the dimensional local mean streamwise ve- 
locity as a function of Reynolds number for Y/d of 0.50, 0.85, and 0.90. Additional 
Reynolds number surveys are given in Figure 11 through Figure 13, which show local 
mean streamwise velocity normalized with respect to bulk velocity, while in Figure 14 
and Figure 15, data are normalized with respect to centerline velocity. In each figure, 
except Figure 11, velocity increases with Reynolds number in the laminar and transi- 
tional regimes with less dependence on Reynolds number in the turbulent region. Such 
increases are largest at Y,d = 0.50. In each regime, at a particular Re, dimensional and 
normalized velocities decrease with increasing Y/d, as expected. In each of the figures, 
data in the transition region from Reynolds numbers of about 2100 to 2500 show some 
e.xperimental scatter. 

2. Longitudinal Turbulence Intensity 

Dimensional longitudinal turbulence intensity data are shown in Figure 16 
through Figure 18. The turbulence intensity data are normalized with respect to 
centerline velocity in Figure 19 through Figure 21 and with respect to local average 
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velocity in Figure 22 through Figure 24. Turbulence intensity data for all three lo- 
cations increase with Reynolds number in the transition region (Re= 1450 to 2500) with 
less dependence on Reynolds number in the turbulent regime when Re becomes greater 
than 2500. Local maxima appear from Re = 2000 to Re= 2400 in each figure which in- 
crease in magnitude with Y'd. In Figure 19 through Figure 21, normalized turbulence 
intensity data show similar qualitative behavior. 

Imposed pulsations at 1 Hz have negligible effect on magnitudes of turbulence 
intensity survey, however with 2 Hz pulsations, intensities with imposed pulsations at 
Y,'d = 0.50 are lower than values without pulsations for Reynolds numbers from 2400 to 
2500 (Figure 22). These variations are not seen at either of the other Y,'d locations 
studied. They are believed to result because of destabilization at 2 Hz resulting in earlier 
onset of fully turbulent flow with less local intensity. 

3. Normalized Phase-Averaged Velocity Magnitudes 

Figure 25 and Figure 26 show variations in peak-to-peak amplitudes of phase- 
averaged velocities normalized with respect to bulk mean velocity (Au/V^rt) as a function 
of Reynolds number. In each of these figures, the data without pulsations is presented 
for comparison with imposed pulsation data. The 1 Hz survey (Figure 25) shows little 
effect of imposed pulsation on Au/V*;*, whereas the 2 Hz survey (Figure 26) shows an 
increase in Au/V^^ with Reynolds number with local maxima apparent near Re =2000 
and Re =2450 for all three Y,d locations. The maxima near Re = 2450 is greatest in 
magnitude. Smaller maxima near Re = 2000 and Re =2450 are also visible in the 1 Hz 
surveys . 

4. Intermittency 

Figure 27 shows turbulence intermittency as a function of Reynolds number for 
each channel Y,'d location of interest at imposed pulsation frequencies ofO, 1, and 2 Hz. 
The intermittency values ( y ) are calculated from instantaneous second derivatives of 
streamwise velocity wuth respect to time. The second time derivative is calculated nu- 
merically using the finite difference method. This value is then compared to a threshold 
value in an on-off temporal identification function, ID, defined such that : 

ID = 1, if > Thr (flow turbulent) 



and : 



ID = 0, otherwise (flow' non-turbulent). 
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Intermitiency is then the percent of the time that the flow is fully turbulent. The inter- 
mittency threshold, Thr, is given by : 

Thr = {105 -(64.66 x (V,„-1.89))} 

This approach produces results as seen in Figure 27. Values for 2 Hz pulsations range 
from 0.5 to 0.8 at low Reynolds numbers (in the laminar region) depending on channel 
location. These values cannot be believed since one expects y to be nearly zero when 
flow is fully laminar. 

In order to provide more reasonable intermittency results, another method for 

computing the turbulence threshold is employed. Since similar intermittency trends are 

evident for each Y/d position, the new threshold is determined at the centerline and then 

taken to be valid for the other positions. To do this, average second time derivatives 
c-U 

of velocity (— ; — ) for the 0 and 2 Hz cases are determined. Figure 28 shows the second 
c t 

time derivative for each case (0 and 2 Hz) multiplied by a factor of 0.01 as a function 
of bulk mean velocity (and consequently Reynolds number). The threshold algorithm 
as described above is also presented in this figure. To determine the threshold factor, 
results for the 0 Hz case shown in Figure 27 are considered to be valid. The value of 
the original threshold algorithm at a given bulk velocity is divided by the average second 
time derivative of velocity to give the new threshold factor, T.F. T. F. values as a 
function of bulk velocity are then shown in Figure 29. A curve is then fit, as shown, to 
express the threshold factor as a function of the bulk mean velocity. This curve is given 
by : 



T.F. = 2.5414VI,,- 13.331V*,, + 17.412 



thus : 



Thr = T.F. X 

ot- 

With this approach, “corrected" intermittency survey results are shown in 
Figure 30 for pulsations at 0 and 2 Hz. In this figure, data for both the 0 and 2 Hz cases 
show intermittency varying from 0 to 1 as Reynolds number increases from 1550 to 
2050. 
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5. Spectral Measurements 

Spectral measurements are presented in Figure 32 through Figure 36 for 
^';d = 0.90 at Reynolds numbers from 1400 to 1700 and pulsation frequencies ofO, 1, 2, 
3, and 4 llz. Figure 31 shows the dynamic signal analyzer (DSA) settings for spectral 
measurement. At Re= 1550, with no imposed pulsations, spectra show "humps" sinrilar 
to ones resulting from spectra of time traces having Dirac delta functions. These corre- 
spond to intermittent turbulent events. As the How is subjected to pulsation, the tran- 
sitional or critical Reynolds number decreases, displaying increasing overall spectral 
energy levels near Re= 1500. 

6. Phase Average Velocity and Turbulence Intensity Traces 

Figure 37 through Figure 76 show phase-averaged velocity and turbulence in- 
tensity traces for Y,'d = 0.85 at 0 and 2 Hz produced as discussed in Appendix A. 
Phase-averaged traces for 2 Hz pulsations (Figure 57 to Figure 76) show sinusoidal 
behavior at low Reynolds numbers which become distorted as Reynolds number in- 
creases. As the bulk mean velocity and local average velocity increase with Re, the im- 
posed pulsations show smaller percent-wise variations on instantaneous velocity traces. 
Increased turbulence levels are more apparent on instantaneous velocity traces and 
phase-averaged traces are increasingly distorted. 
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V. SUMMARY AND CONCLUSIONS 



This study examines laminar, transitional, and turbulent flows in a 40 to 1 aspect 
ratio straight channel at Reynolds numbers from 1 100 to 3400. The study also investi- 
gates the influences of pulsations imposed on the flows at frequencies of 1, 2, 3, and 4 
Hz. Data are presented for Y/d locations of 0.50 (centerline), 0.85, and 0.90. Stokes 
numbers are 4. OS and 5.79, and Strouhal numbers vary from 0.039 to 0.122. 

Local average velocity surveys, both dimensional and normalized, are largely unaf- 
fected by imposed pulsations. Longitudinal turbulence intensity levels for the three lo- 
cations show increases with Reynolds number from Re= 1100 to Re = 2200, and then 
show decreases with Reynolds number from 2200 to 2500. Afterwards, smaller vari- 
ations with Re are apparent as the flow becomes fully turbulent at Re >2500. Longi- 
tudinal turbulence intensity values which are higher at the channel centerline than near 
the wall evidence a center mode of transition for Re from 1200 to 1500 (Zang and Krist 
(Ref 28]). Maximum turbulence intensity values observed near Re= 2200 increase with 
Y/d, as does the turbulence intensity in the fully turbulent region. 

For Reynolds numbers from 2350 to 2450, imposed pulsations at a Stokes number 
of 5.79 (2 Hz) result in decreased normalized turbulence intensity values compared to 
values seen at 0 and 1 Hz. This is believed to occur because pulsations destabilize the 
flow, an effect not seen for 1 Hz pulsations, perhaps because of their lower amplitudes. 

Spectral measurements show overall energy level increases to occur between Re of 
1500 and 1550, depending on pulsation frequency. Spectra show evidence of intermit- 
tency at Re of 1550 for imposed pulsations at 0, 1, 2, 3, and 4 Hz. Intermittency values 
for no imposed pulsations and for 2 Hz imposed pulsations increase from 0 to 1 as Re 
increases from 1450 to 2100. 

Phase-averaged velocity traces show sinusoidal behavior over low’ Re below 1750. 
As the Reynolds number increases, sinusoidal traces become more distorted at Re of 
1850 to 2200 in the transition regime, and at higher turbulent Re W'hen turbulence in- 
tensity levels increase. Amplitudes (peak-to-peak) of phase-averaged traces show im- 
portant variations with Re. With 2 Hz imposed pulsations, local maxima are evident 
near Re of 2000, 2200, and 2450. 
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APPENDIX A. SOFTWARE DIRECTORY 



DEAN 15 : Determines channel mass flow rate, bulk velocity and Reynolds number 
from orifice plate pressure drop. The program requires ambient temperature and pres- 
sure inputs which can be entered manually or via a data acquisition system connected 
to thermocouples and/or pressure transducers. Once ambient conditions are established, 
air density is calculated and orifice pressure drop is entered manually. Originally de- 
signed for curved channel applications, this version has been specifically modified for 
straight channel use. 

HWCAL : Determines the constants for the King's Law calibration of the hot-wire. 
The program also provides a polynomial fit of the calibration data. The version used is 
written for the IB.M 310 mainframe computer. 

HOTAVIREPAV : Reads the data stored in the A/D buffer of the high speed data 

acquisition system and stores the information on micro diskettes. Manual inputs are : 
triggering frequency, hot-wire DC voltage (ungained), oscillation frequency (flow 
blockage), bulk velocity, Reynolds number, and date and time of run. 

DATAP : Calculates instantaneous and phase-averaged velocities. Initially a look 
up table is created. Effective velocities are calculated from the effective voltage values 
and stored for follow-on calculations. The hot-wire calibration constants obtained from 
HWCAL, and the amplifier gain are incorporated into these calculations. The velocity 
calibration is given by : 



= k{Ei,- Eir^ 

where k is the proportionality constant, E.^ is the effective voltage, and Eo is the refer- 
ence voltage at no flow. N is a constant (0.45) for moderate Reynolds numbers. Once 
the look up table is created, the program reads the instantaneous voltage values from the 
data file and converts them to instantaneous velocities. A plot of the instantaneous ve- 
locity versus time can be generated for any of 580 cycles (the first 60 cycles are discarded 
to allow for flow stabilization as the unsteady device is started). 

Next the program phase averages the 580 cycles, and velocity versus phase angle (of 
the flow blockage) plots arc available. Two plots are available from the averaged 
values : u versus phase angle and /u versus phase angle, u is the phase averaged 
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velocity, yj u'^ is the phase averaged root mean squared (rms) velocity, and u is the av- 
erage velocity. In the case where there is no imposed unsteadiness, the phase averaged 
velocity, u , is equal to the time averaged velocity, u. 

INTRMTCY : Takes data recorded by HOTWl REPAY and computes the inter- 

mittency. It uses the second time derivative of the instantaneous velocity, which is cal- 
culated numerically using the finite difference method. This value is then compared to 
a threshold value in an on-off temporal identification function, ID, defined such that ; 

ID = 1, if 4^ > Thr, and 
cv 

ID = 0, otherwise. 

The threshold factor, Thr, is determined for different Reynolds numbers using the 
correlation : 



Thr = {105 -(64.66 x (V,,* - 1 .89))} 

If the Thr value falls less than 20 then its value is reset to 20. The values of ID are then 
phase averaged and time averaged for the traces and mean intermittency value, y . This 
approach produces results which are reasonable for all but the 2 Hz imposed pulsation 
case. 

In order to provide more reasonable results for intermittency, the method for com- 
puting the turbulence threshold is modified. Since the intermittency surveys show a 
similar trend for each Y/d position, the modification approach for the threshold is done 
at the centerline and taken to be valid for the other positions. The new thres hold a ko- 
rithm requires investigation of the average second time derivative of velocity ( ., ^ ) for 

C/ X 

the 0 and 2 Hz case. The goal is to obtain a "threshold factor" to multiply the second 
time derivative since the intermittency is a function of the magnitude of the second time 
derivative as well as the bulk Reynolds number. 

The approach for determining the threshold factor assumes the results for the 0 Hz 
case to be valid. The value of the original threshold algorithm at a given bulk velocity 
is divided by the average second time derivative of velocity. This is the threshold factor, 
T.F. A curve is then fit to express the threshold factor as a function of the bulk mean 
velocity. Now : 



Thr 



T.F. X 
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where : 



T.F. = 2.5414V^*- 13.331V*,, + 17.412 



Since the threshold modification approach produces reasonable results for the 0 Hz case, 
which had previously been believable, it is assumed that the modification of Thr will 
have no significant effect on the 1 Hz case. 
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APPENDIX B. 
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Figure 1. Straight Channel Schematic 
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Figure 3. Phase Angle and Mean Velocity Variation nith Vane Position 
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SCHEMATIC OF FLOW MEASUREMENT EQUIPMENT 




Figure 4. Flow Measurement Equipment 
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Figure 5. Bulk Velocity vs Orifice Differential Pressure 
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Figure 6. Channel Pressure Distribution 
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Figure 7. Inlet Turbulence Intensity Conditions 
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Figure 10. Vavg vs Re (Y/d = 0.90) 
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Figure 1 1. Vavg/Vbik vs Re (Y/d = 0.50) 
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Figure 12. Vavg/Vbik vs Re (Y/d= 0.85) 
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Figure 13. Vavg/Vbik vs Re (Y/d= 0.90) 
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Figure 14. Vavg/Vcl vs Re (Y/d = 0.85) 
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Figure 15. Vavg/Vcl vs Re (Y/d = 0.90) 
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Figure 16. Vrms vs Re (Y/d = 0.50) 
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Figure 17. Vrms vs Re (Y/d= 0.85) 
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Figure 18. Vrms vs Re (Y/d= 0.90) 
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Figure 19. Vrms/Vcl vs Re (Y/d= 0.50) 
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Figure 20. Vrms/Vcl vs Re (Y/d= 0.85) 
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Figure 21. Vniis/Vcl vs Re (Y/d= 0.90) 
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Figure 22. Vnus/Vavg vs Re (Y/d = 0.50) 
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Figure 23. Vriiis/Vavg vs Re (Y/d= 0.85) 
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Figure 24, Vrm.s/Vavg vs Re (Y/d = 0.90) 
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Figure 26 . Delta Uhat/Vbik vs Re (2 Hz) 
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Figure 27. Intermittency vs Re (1 Hz, 2 Hz) 
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Figure 28. Threshold and Second Time Derivatives vs VbIk (Y/d = 0.50) 
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Figure 30. Corrected Inleriiiittency vs Re (2 Hz) 
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Figure 31. Typical Dynamic Signal Analyzer Settings 
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Figure 32. Spectral Frequency Response (0 Hz) 
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Figure 33. Spectral Frequency Response (1 Hz) 
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Figure 34. Spectral Frequency Response (2 Hz) 
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Figure 35. Spectral Frequency Response (3 Hz) 
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Figure 36. Spectral Frequency Response (4 Hz) 
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Figure 37. Average and RMS Velocity Traces (Y/d = 0.85, Re = 1100, 0 Hz) 
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Figure 38. Average and RMS Velocity Traces (Y/d = 0.85, Re= 1250, 0 Hz) 
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Figure 39. Average and RMS Velocity Traces (Y/d= 0.85, Re= 1350, 0 Hz) 
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Figure 40. Average and RMS Velocity Traces (Y/d= 0.85, Re= 1440, 0 Hz) 
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Figure 41. Average and RMS Velocity Traces (Y/d = 0.85, Re= 1550, 0 Hz) 
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Figure 42. Average and RMS Velocity Traces (V/d= 0.85, Re= 1710, 0 Hz) 
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Figure 43. Average and RMS Velocity Traces (Y/d= 0,85, Re= 1850, 0 Hz) 
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Figure 44. Average and RMS Velocity Traces (Y/d = 0.85, Re= 1920, 0 Hz) 
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Figure 45. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2000, 0 Hz) 
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Figure 46. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2040, 0 Hz) 
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Figure 47. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2090, U Hz) 
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Figure 49. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2300, 0 Hz) 
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PHASE mZLZ (DEC.> 



G5_2350_O 10391 . I 156 

RVG. VEL 2.-A53mx^3 RMS VEL 

OSC. FREO. 0Hz STR. # 
BULK VEL. 2.8 1m-'£ REY. ♦ 



^ ^ 2 1 m 3 
0 

2350 




85_2350_0 10391 . 11 5B 

fS\J6, 1;EL 2.459m/5 RMS UEL .4421n/s 
OSC. FREO. 0Hz STR. « 0 

8ULK OEL. 2.81n/5 REY. t 2350 

Figure 50. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2350, 0 Hz) 
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05 _ 2 ^ 00_0 1 12790. 1951 

RVG. VEL 3.016m^s RMS VEL .501Sm/3 
CSC. FREO. 0H2 STR. ^ 0 

BULK VEL. 2.072m.-s REY. # 2400 




Pfiflse RNOLC <k:c.> 



95_24O0_G 1 12790. 1951 

RVG. VEL 3.016m/3 RMS VEL .501Sm/3 
OSC . FREO. OMz STR. 0 

BULK VEL. 2.872m-'"s REY. # 2400 




85_2400_0 I 12790.1951 
A^G. VEL 3.0lBn/6 RMS VEL .50l9ri/5 
OSC. FREO. 0Hz STR. t 0 
bulk VEL. 2.872n/5 REV. t 2400 



Figure 51. Average and RMS Velocity Traces (V/d= 0.85, Re = 2400, 0 Hz) 
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e ^ 5 0_0 1089] .2244 

nVG. VEIL 2.735m/s RMS VEIL .45B9m-^c 
OSC . FREO. 0Hz STR. # 0 

BULK' VEL . 2. 93 5m REV. +► 2 45 0 




8 5_2 4 3 0_0 10591 .2244 

RVG. VEL 2.7 3 5rr..^» RMS VEL . 4 58 9m^s 
OSC. FPEO. 0Hz STR. 44 0 

BULK VEL. 2.935m-'s PE'V . 44 2450 




PHASE mCLC (DEO.) 



85_2450_0 10891 .2244 

AUG. UEL 2.735r./s RMS UEL .4589n/s 
OSC. FREO. 0Hz STR. # 0 

BULK UEL. 2.935n/5 REY. t 2450 



Figure 52. Average and RMS Velocity Traces (Y/d= 0.85, Re = 2450, 0 Hz) 
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. ^54 3m-^s 



B3_2300_O 120130.221 

f=<VG. VEL 2.9B3m>^s RMS VEL 

05C. FREO. GHz STR . # 0 



PULK VEL. 2, 99 2m.'- s REV. -♦+ 2 500 




FWRSE RJCLX (DEG. ) 



B 3 _2 5 00_0 12 0 13 0.2213 

HVG. VEL 2 . S 0 3r.-,."s RMS VEL .4 54 3n-,.'s 
DSC. FREO. 0Hz STP. # 0 

BULk VEL. 2.992m-'s PE i . -♦+ 2500 




PVJfliC RIGLC (DCO.i 



c5_Z5O0_0 120190. 2Z15 

Ai;6, VEL 2.8B5n/5 RMS VEL .4543n/s 
nSC. FREQ. 0Hr STR. t 0 
BULK VEL. 2.992n/s REY. « 250G 



Figure 53. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2500, 0 Hz) 
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RR_?B(?l(?(_C?» J?C(E9Pi.?25g 

ncn. FRFo, 0H7 ftr. « 0 

PULF VEL : rev., t 7R?n 



Figure 54. Average and RMS Velocity Traces (V/d= 0.85, Re = 2800, 0 Hz) 
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B3_3 1 00_0 1043] .2133 

RVG. VEL 3.327mv-s RMS V'El_ . 4 7 1 1 rr,--s 
OSC . FREO. 0Hz 5TR. # 0 

BULK VEL- 3.7 1 Im/s REY . 4f 3 100 




85_3 1 0O_0 10431.2135 

RVG. VEL 3.327m^s RMS VEL . 4 7 1 1 m^s 
OSC. rPEO. 0 Hz STR. # O 
BULK VEL. 3.71 Im^s RE>' . # 3 100 




PWRSC R'4GLC (DCG.) 



85_3100_0 10491 .2155 

AVG. VEL 3.327n/5 RMS VEL .47Mn/5 
OSC. FREO. 0H2 STR. t 0 
BULK VEL. 3.7lln/5 REY. t 3100 



Figure 55. Average and RMS Velocity Traces (Y/d= 0.85, Re = 3100, 0 Hz) 
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B5_3 5 5 5_0 3 3 I I B O . 05 5 2 

PVC. v’CL- 3-324m-'s RMS V’EIU -4 55 5n-»''s 
DSC- r^EO. 0Hz STR . 0 

PULK VEL. PEIv, 4f 2 5P5 




PMPSC mcLC (DCG. ' 



RVC 



VEL 3 
09 C . 



5 =; p «5 rn 

,334 n*i-* s 

F K El O . 



3 3 I I B 0 0S 3 
RMS v’ZL .4 5: 
0 »4 2 9 ^ ^ 0 



PL'L^ V/0L-, 



PE »' . 




t- 



‘ V« */*t »efc> mV \et 2 11- ^ zbi kb6 S'i'^ 39's' 

FVfOSC «=fVGLC ^DCC. ' 



J?9_Zbi9 1 t 1 ] 9i? . l?bbZ 

^••u iiCi 7 7 7, ^^/C’ Rf1'“ 

ObC. FREQ. 0Hr ETR. 5 C 
SL’LK I'EL. 4.t?49"^./£ REV. 4 



4555n ' £, 



Figure 56. Average and RMS Velocity Traces (Y/d = 0.85, Re = 3400, 0 Hz) 
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B5 1 L 30 2 



} : 1 0SO . 17]] 



RVC. VE:i_ .7592n^--« RflS 

OSC . FRlO. 2 Hz STP. # 
PUL^^< VEL. I . ?09n>-'s‘ REy 



00 2 5 05 K. 
- I 2 I 5 
^ 115 0 



5 




tok- ^ 

Kkt I * I .... t ^ 

38 7 k i&b i 6 b sit eSk ktb see 



FVPSE FTiCtE fECG.J 



® * 2 13 0 

R v'C . v'El_ . 7 5 9 2 ft’t'- G. 

Or-C . F E O . Z Hz 
PUL^^- . 1 . 509-' 



2 2 I 05 0 .17 2 2 
RMS v'El_ . 002 5 05 n*.- s 
r • T F . # . I Z I 5 






115 0 




Fwnsc WCLC fDCO.» 



fib_l 1 1 I 1 e9i?. 1 7M 

MVt. v£l .7bS2r.,^£ Rn£ f£L .Cl?250£fn.^s 
OSC. FREO. 2hz STR. t .I2i9 
BULK VEL. ].3e5r,.'£ REV. « MbO 



Figure 57. Average and RMS Velocity Traces (Y/d = 0.85, Re= 1100, 2 Hz) 



76 



9 S_ I 2 5 0_2 122090. 15-19 

RVG. VEL .9533m^s RMS VEL .003596m/3 
CSC. TREO. 2Hz STR. # .3071 

BULK VEL. REY. # 1290 



or 




05 1250 2 I2299Q. 15-19 



RVG. VEL .9533m/s RMS VEL 
CSC. FREQ. 2H2 STR. ^ 
BULK V/EL. R E i' . 



. 003596m/3 
. 3071 
♦ I c 9 0 




85_1250_2 122890.1549 

AUG. UEL .9513m/s RMS UEL .003596n/s 
OSC. FREO. 2Hz STR. t .1071 
BULK UEL. 1.49r./s REY. t 1250 



Figure 58. Average and RMS Velocity Traces (Y/d = 0.85, Re= 1250, 2 Hz) 



77 



122 - 190 . 102 ? 



05_l 350_2 

RVG. VEL 1.05?n./3 RMS VSL .0l?e7m^s 
03C. FRSCj. 2H2! sir. It .0303? 
BULK V'EL. l.622rr\.-'5 RSY . # 13^0 

f- 

s> 




05_1350_2 122-190. 1027 

RVG. VEIL 1.057m-'3 RMS VEL .0l7e7m.^s 
05C . FREO. £H2 STR . It .0385? 
BULK VEL. l.622nu^s RET. It 13!5G 




85_1350_2 122490.1827 

AUb. UEL 1.057r>/a RMS UEL .01787n/5 
OSC. FREQ. 2Hz STR. t .09337 
bulk UEL. 1.622n/s REY. I 1350 



Figure 59. Average and RMS Velocity Traces (Y/d= 0.85, Re= 1350, 2 Hz) 
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Figure 60. 



B5__l 5G0__2 2 1 1390.3352 

RVC . VEL .95n..'« Rt IS VEL .Q2\27n-xys, 
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CSC. FREQ. 2 hz ETR. t .C9293 
Bulk V£l. 1.7l7m/s REV. > I see 



Average and RMS Velocity Traces (Y/d = 0.85, Re= 1440, 2 Hz) 
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G5_l 5 5 0_S 120090.2205 

RVG. VEIL 1.16001/3 RMS VEL .0506?m/3 

05C. FREQ. 2Hz 5TR. .0658 

BULK VEL. l.eSm/s REY . ♦ 1550 

m 




©5_l 5 5 0_2 12 00 90.2205 

RVG. VEL 3.16001/5 RMS VEL .05067m/3 
05C . FREO . 2Hr STR. ♦ .0658 

BULK VEL-. l.eSm/s REY. ♦ 1550 




85_1550_2 120890.2205 

AUG. VEL l.lB8m/5 Rn5 VEL .05067n/s 
05C. FREQ. 2Hz STR. t .0858 
BULK VEL. 1.86n/5 REY. t 1550 

Figure 61. Average and RMS Velocity Traces (Y/d = 0.85, Re= 1550, 2 Hz) 
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05 l?50 2 



1 1 1590.0109 



RVG. VEL 1.239m/3 RMS VEL -115<4m/3 
05C . FREO. EHz STR. # .0?679 

BULK VEL. 2.07em-'-s REY . ^ 1750 




0 5_l 7 5 0_3 
RVG. VEL 1 . 239m/3 
05C . FREO. 2 Hz 
BULK VEL. 2.07Bm. 



1 1 1590.0109 
RMS VEL . ] lS*^m/3 
STR. ^ .07679 

REY. 1750 




85_1750_2 n 1590.01C9 
A^6. ^El 1.239rn/s RhS VEL .1154n/s 
OSC. FREQ. 2Hz STR. t .07679 
BULK VEL. 2.078m/s REY. t 1750 

Figure 62. Average and RMS Velocity Traces (Y/d= 0.85, Re = 1710, 2 Hz) 



81 



G5_10 5O_2 12(3090.1013 

RVG. VEL 2.-195m/3 RMS VEL. . 2 1 9 1 m/a 
05C. FREcj. 2Nz STR. m .0?217 
BULK VEL. 2.21 Im^s REY. # 1050 



CM 

ir> 




PHRse W4CIX <k:c.> 



65_1950_2 120090. 1013 

RVG. VEL 2.-^95m/3 RHS VEL .2191m/3 
05C. FREO. 2Mz STR. n .0721? 
BULK VEL. 2.2llnt.'-s REY. n 1050 




65_1 850^2 120890.1813 

AUb. UEL 2 . 495 p »/5 RMS UEL .2l91n/s 
ObC. FREO. 2Hz STR. t .07217 
BULK UEL. 2.2Mn/5 REV. t 1850 

Figure 63. Average and RMS Velocity Traces (Y/d= 0.85, Re= 1850, 2 Hz) 
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05__2005_2 I 1 1 090 .210? 

RVG. VEL 1 . 653m.'s RMS VEU -303Sm.-3 
CSC. FREO. EH:: 5TR. # .0637 1 

BULK VEL. 2.2S9rT>."3 PET. 4t 2005 




©5 2005 2 111090.2107 



R v’G . VEL 1 . 6 5 3 m 3 RM E L . 3 O 3 & m 3 

05C . FREO. 2Hr 5TR. ^ .06971 

BULK VEL. 2. 2 6 9m. "3 PEY. 2005 




85^2005.2 111090.2107 

AVG . UEL 1.653n/s RMS UEL .3059n/s 
05C. FREO. 2Hz SIR. f .06971 
BULK VEL. 2.2e9r./s REY. t 2005 

Figure 64. Average and RMS Velocity Traces (Y/d= 0.85, Re= 1920, 2 Hz) 
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5 I eO_2 J 1 1 I 50 . ce } 2 
RVC. VCL RMS VEU . 3325n-.^'e 
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Figure 65. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2000, 2 Hz) 
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0 5__2 3 5 0_2 1 1 1 I 9 O . 02 5 S 

BvG. VEL. l.eGlm.'3 RMS VEL. .3564^^3 
CSC. FREQ. 2Hi: STR. .065^3 

BULK VEL. 2. 4 37m ■" -5 REY. 4 2 150 







e 5_2 3 5 0__2 1 1 1 I 9 0 . 02 5 G 

RVG. VEL l.eeim.'r RMS VEL .3564m.'r 
OSC. FREO. 2Hz STR. .06545 

BULK VEL. 2.4 37m.''s REY. 2 150 




Figure 66. 



e5_2i50_; 111190.0:56 

AUG. UEL 1.86lr>/s RMS VEL .3564n/s 
OSC. FREO. 2Hz STR. # .06549 

BVLF VEL. Z.457n/s REY. # 2150 



Average and RMS Velocity Traces (Y/d = 0.85, Re = 2040, 2 Hz) 
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112090.0212 



G 5 2200 20 

RVG. VEL_ 2.07 1m/3 RMS VEL. .-4224m/3 
03C . TREO. STR. ft .06066 

BULK VEL . Z.63lm-"s REY . ft 22 00 




05_22OO_2O I 12090. 02 12 
RVG. VEU 2.07 1m/3 RnS VEU .^22-^m/3 
CSC. FREO. 2Hz STR. ff .06066 
BULK VEL. 2.63lm.-'s REY. ft 22 00 




85.2200,20 112090.0212 

AU6. UEL 2.071n/s RMS UEL .4224m/5 
OSC. FREO. 2Hz STR. t .0B0BB 
BULK MEL. 2.B3ln/5 REY. t 2200 



Figure 67. Average and RMS Velocity Traces (Y/d= 0.85, Re = 2090, 2 Hz) 
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123090. 1005 

RVG. VCL 2. 533m/3 Rf^S VEl_ -^530m/3 
CSC. FREQ. 2Hz STR. ^ .0605^ 

BULK VEL. 2.636m--s REY . # 2200 




PHRSE RNCLE <DEG.) 



G5_2200_2 123090. 1905 

R VG . V E L 2 . 5 3 3 m v' r Rl"1 5 VEL .^530nn^3 
CSC. FREQ. 2Hz STR. # .0605-=^ 

BULK VEL. 2.636m.-s REY. 2200 




Figure 
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05_2200_2 123090.1605 

AVG. VEL 2.533n/s RMS UEL .4S3Bn/s 
OSC. FREO. 2Hz STR. # .06054 

BULK UEL. 2.636n/5 REY. I 2200 



Average and RMS Velocity Traces (Y/d= 0.85, Re = 2200, 2 Hz) 
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0 5_2 3 00_2 1 I 2G90 . 22 I 1 

f=IVG . VEL 2.72m/3 RMS VEU .5l-»3m^3: 
OSC. FRECr. £ Hz 5TR. It .05801 
BULK VEL. S.rSlm.'^s REY. tt 2300 




PH«3C W4CLE (DGC.) 



0 5_2 3O0_2 
BVG . VEL 2 . 72m/3 
OSC. FREO. 2Hz 
BULK VEL . 2 . *75 I m 



1 1 2G9C . 22 I 1 
RMS VEL .5l^3m^» 

5TR . It ,0560 1 

.-s REY. tt 2300 




85.2300.2 112690.2211 

htJG. UEL 2.72n/s RMS VEL .5143fn/s 
OSC. FREO. -2Hz SIR. * .05801 

BULK UEL. 2.751n/s REY, I 2300 



Figure 69. Average and RMS Velocity Traces (Y/d= 0.85, Re= 2300, 2 Hz) 
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O5_2350_2 10391 . I G 1 9 

RVG. VEl_ 2.-^?2m/3 RMS VEL 

CSC. FREQ. 2Hz STR. ft .0567S 
BULK VEL. a.eim^s REY . -ft 2350 




PHASE ANv:lE (DEG. ) 



9 5 _2 3 5 0_2 10 3 9 1 . I G I 9 

R VG. VEL 2.-A72m/s RMS VEL -ft -A 60^/3 
CSC. FREO. 2H2: STR. ft .05673 

BULK VEL. 2.8 1m^'s REY. ♦ 2350 




85_2350_2 10391 . IB 19 

AUG. VEL 2.472n/5 RMS VEL .4488n/s 
OSC. FREO. 2Hz STR. # .05679 

BULK VEL. 2.8ln/s REY. t 2350 

Figure 70. Average and RMS Velocity Traces (V/d = 0.85, Re = 2350, 2 Hz) 
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0 5_2^Q0_2 1 12090. 170-1 

RVG. VEL. 2.670m/3 RMS VCL_ .-16©4m/3 
05C. FREO. 2Hz STR. n .05549 
BULK VEL. 2.0?6m-^s REY , # 24 00 




PMRSC R4GLC (ICC. > 



0 5_2 4O0__2 1 120 90. 17 04 

RVG. VEL 2.670m/3 RMS VEL -46e4m/s 
03C. FREO. 2Hz STR. It .05549 
BULK VEL. 2.876m.^2: REY. It 2400 




85^2400^2 112890. 1704 

AUG. UEL 2.B78n/s RMS UEL .4G84n/s 
05C. FREO. 2Hz STR. t .05549 
BULK UEL. 2.878^/5 REY. t 2400 



Figure 7L Average and RMS Velocity Traces (Y/d = 0.85, Re = 2400, 2 Hz) 
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B ^ 3 0_2 I eJB 9 3 .2352 
RVC. VEIL 2.75lm/s RMS VEL .4^3 1m--s 
OSC. FREO. 2Hz STR . # .05438 

BULK VEL. 2.935m/s REV. # 2450 




PMPSC RJCLE (DCG.) 



B5_2 4 3 0_2 
RVC . VEL 2-751 m-''s 
OSC, FREQ. 2Hz 



10831 .2332 
RMS VEL .4431m- 
STR. # .05438 



BULK VEL. 2.935m-^s 



PEY . 4f 



2 45 0 




85_2450_2 10891 .2352 

AUG. UEL 2.75ln/s RMS UEL .4431tn/s 
OSC. FREQ. 2Hz STR. t .05438 
BULK UEL. 2.935rn/5 REY. t 2450 



Figure 72. Average and RMS Velocity Traces (Y/d= 0.85, Re = 2450, 2 Hz) 
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63_2300_2 120190.2316 

RVG. VEIL 2.809m/« RMS VEIL .4502m^s 
CSC. FREO. 2Hz STR . # .05334 

BULK VEL. 2.992m^s REV. # 2500 




n>iCLC CDCG.) 



6 5_2 5 00_2 
RVG. VEL 2.899 s 
CSC. FREQ. 2Hz 
BULK VEL . 2.992m 



120190.2316 
RMS VEL .45G2rr,^'s 
STR. # .05334 

-s REY. 4f 2 5 00 




R-JOLC CDCO.> 



85^2500.2 120190. 251B 

AV6. VEL 2.889n/5 RMS VEL .4502n/s 
OSC. FREO. 2Hz STR. « .05334 

BULK VEL. 2.992n/s REY. t 2500 



Figure 73. Average and RMS Velocity Traces (Y/d=0.85, Re = 2500, 2 Hz) 
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83_2000_S 120?80.00O3 

RVG. VEIL 2.B21m^s RMB VEIL .3983m^s 
CSC. FREO. 2Hz STR . # .04782 

BULK VEL. 3.35 Im^s PEY . # 2000 




FV»«SE mCLC (DEC.) 



B5_2 8 00__2 120790.00 03 

RVG. VEL Rr^S VEL .398 3m.'s 

05C. FREQ. 2Hz STR. # .04762 

BULK VEL. 3.35 1m.^s REV. # 2900 




PMncC PTiCtC (DEC. ) 



B5_2800_2 120790.0003 

^K}G. UEL 2.821n/s RMS UEL .39B3rn/s 
OSC. FREQ. 2Hz STR. t .04762 
BULK UEL. 3.351n/5 REY. t 2B00 



Figure 74. Average and RMS Velocity Traces (Y/d = 0.85, Re = 2800, 2 Hz) 
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Figure 15 . 



B5_3 1 00_2 1 3 } .2 304 

RVC. VE:1_ 3.332m-^s RMS VEL. .4533m/s 

05C. FREQ. 2Hz 5TP. # .043 

BULK VE1_. 3 . 7» 1 Im/® REY . 4f 3 100 




mnsc mcLC (Dtc. ) 



B5_3 100_2 1043 1.2304 

BVG. VEL 2.3 3 2m^-s RMS VEE .4 53 3rr,.-s 
OSC. FREO. 2Hz STR . # .043 

BULK VEL. S.'Pllrr,-^® REN'. 4f 3 100 




PWflSC mCLC <DCG.) 



85_5100_Z 10491 .Z304 

AUG. UEL 3.352fn/s RMS UEL .4553fn/s 
OSC. FREO. 2Hz STR. t .043 
BULK UEL. 3.71 ln/5 REY. t 3100 



Average and RMS Velocity Traces (Y/d = 0.85, Re= 3100, 2 Hz) 
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05 _ 3505_2 1 1 1 190.0907 

RVG. VEL 3.26m/3 RMS VEL .45-l9m^:s 
03C. FREO. 2HZ STR. ^ . 039-13 

BULK VEL. -1.0-10m.'"s REY, ^ 3385 



r- 

OP 




05_3505_2 1 1 1 190.0007 

RVG. VEL 3.26m/3 RMS VEL .H5-»9m^3 
05C. FREO. 2Hz STR. # .039*13 

BULK VEL. -1.0*iem--s REY. # 3385 




85_3585.2 I 1 1190.0807 
WG. l^EL 3.2Bn/5 RMS VEL .4549n/s 
OSC. FREQ. 2H2 STR. t .03943 
BULK VEL. 4.048n/s REY. t 3585 

Figure 76. Average and RMS Velocity Traces (Y/d= 0.85, Re = 3400, 2 Hz) 
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